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1 Introduction 

In the vast and diverse field of quantum magnetism, the quantum properties of large 
magnetic molecules have enjoyed a strong and motivated research activity for more 
than a decade. Physicists and chemists, theoreticians and experimentalists, engi- 
neers and philosophers, all would find at least one good reason to be interested in 
these systems. We shall refer here to "quantum nanomagnets" as the broad family 
of molecules containing a core of magnetic transition-metal ions, which interact by 
strong superexchange and which possess magnetic anisotropy due to crystal field 
effects lfl3l . At sufficiently low temperatures (typically ~ 10 K), i.e. much lower 
than the typical intramolecular exchange interaction energy, the whole molecule ef- 
fectively behaves as a nanometer-sized magnet. The total ground state spin value 
can be rather large, S ~ 10, and is often called "giant spin". The best studied exam- 
ples are Mni2 |22|. Fes [50| and Mn4 |3]. The magnetic core of each molecule is 
stabilized by organic ligands, and the molecules are typically bound to each other by 
van der Waals forces to form electrically insulating crystals (Fig. [TJ. With just a few 
exceptions |49l[T6l . the magnetic interaction between different molecules is only of 
dipolar origin, thus orders of magnitude weaker than the intramolecular exchange. 
Single-ion anisotropies, exchange interactions, point symmetry and crystal struc- 
ture all contribute in a complicated way to the total magnetic anisotropy of the giant 
spin. The first great success of molecular magnets consisted in the demonstration of 
magnetic bistability and hysteresis at the molecular level ll35l . i.e. not arising from 
long-range interactions but only from local ones. Molecules possessing this sort of 
bistability were named "Single Molecule Magnets" (SMMs). One could then think 
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Fig. 1 Crystalline and molecular structure of a prototypical quantum nanomagnet, Mn^-ac 1221 
[36). The core contains 8 Mn 3 + (s=2) and 4 Mn 4 + (S=3/2) ions, which interact by superexchange 
to form a ferrimagnetic ground state with total spin S=10. The molecules crystallize in a tetragonal 
structure. 



of 2D arrays of such molecules QUI as the ultimate magnetic recording medium 
fl8l . with information density ~ 10 12 bit/cm 2 . For this purpose, the main research 
goal is to achieve the highest possible single-molecule anisotropy, to allow magnetic 
bistability up to - ideally! - room temperature. 

However, a proper description of bistability in nanomagnets would be incom- 
plete without quantum mechanics. Since the anisotropy barrier for reorientation of 
the giant spin is large but not infinite, there is in principle a non-vanishing probabil- 
ity for the spin to invert its direction by quantum-mechanical tunneling through the 
barrier [7|. That is, the magnetic memory would delete itself because of quantum 
mechanics! The tunneling probability can be estimated with the knowledge of the 
anisotropy parameters of the giant spin, obtained e.g. by electron paramagnetic res- 
onance or neutron scattering. The resulting tunneling rate turns out to be extremely 
sensitive to the external magnetic fields applied to the molecule. In Mn^ for in- 
stance, one would naively expect the tunneling rate to become astronomically long 
if a stray field of just 10~ 9 T is applied along the anisotropy axis. The experimental 
observation of quantum tunneling of magnetization in Mni2 BT1 [lOl [T31 [121 repre- 
sented a major breakthrough, but in some sense a puzzling one. 

The puzzle was solved by carefully considering the role played by the coupling 
between giant spin and surrounding nuclei, which are always present in the ligands 
I3 C, 35 C1, . . . ) or in the magnetic ions themselves ( 55 Mn, 57 Fe, 53 Cr, . . . ). The 
dynamics of the nuclear spins generates a fluctuating magnetic field on the giant 
spin, thereby sweeping its energy levels through the tunneling resonance and yield- 
ing a finite tunneling probability ll32ll . Once the giant spin is allowed to have (quan- 
tum) fluctuations of its own, it exerts a back-action on the nuclei and essentially de- 
termines their dynamics. Thus, the system of giant spin + nuclei must be analyzed as 
a whole, and the resulting theoretical description is now known as "theory of the spin 
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bath" [ 34-| . It predicts, among other things, a square-root law for the relaxation of 
magnetization at very low-T (331, an d a dependence of the tunneling rate on the iso- 
topic composition of the sample. These predictions have been promptly verified by 
measuring the electron spin relaxation in Fes crystals 1 46 48 1 , where it is possible to 
obtain samples that have stronger (by 57 Fe enrichment) or weaker (by replacing *H 
with 2 H) hyperfine couplings as compared to those with natural isotopic abundance, 
while leaving the giant spin unchanged. The picture has been completed by experi- 
ments looking at the "other side of the coin", i.e. studying directly the dynamics of 
the nuclear spins by NMR experiments 03[n][IIl|25]|27l|4]], and finding that such 
dynamics is indeed profoundly entangled with the quantum fluctuations of the giant 
spin. At the present stage, much of the efforts in the field are directed towards using 
quantum nanomagnets as spin qubits for quantum information purposes, by pushing 
the giant spins into a regime where the tunneling can be made coherent ll29ll . 

The purpose of this mini-review is to give an accessible and somewhat pedagog- 
ical introduction to the crucial aspects of the coupled system "quantum nanomagnet 
+ nuclear spins", what makes it special, what has been understood, and what requires 
further attention. It will be shown that much of the common knowledge on nuclear 
and electron spin dynamics is entirely inappropriate to describe this system. Impor- 
tantly, the same often applies to other quantum degrees of freedom in mesoscopic 
physics (SQUIDs, quantum dots, . . . ). The discussion given here on quantum nano- 
magnets should thus be taken as the "worked-out example" of a problem of much 
broader interest. The reason for choosing this specific example lies in the wonder- 
ful property of quantum nanomagnets to combine mesoscopic size and fascinating 
physics, with the cleanliness and reproducibility of a product of synthetic chemistry. 



2 Theoretical framework 

At temperatures such that the internal magnetic excitations can be neglected, i.e. 
when the entire molecule behaves as a giant spin of value S, the minimal effec- 
tive spin Hamiltonian that describes the quantum nanomagnet coupled to a bath of 
nuclear spins {4} is: 

J? = -DS 2 Z +E(S 2 - Si) - £MbS ■ B - £l,A,S, (1 ) 

k 

where z is the easy axis of magnetization. The first term of this Hamiltonian, —DS 2 , 
gives rise to an energy levels structure as shown in Fig. |2a), with a parabolic 
anisotropy barrier separating pairs of degenerate states. Were this the only term 
in the Hamiltonian (HJ, its eigenstates would be the eigenstates of S z , i.e. the \m) 
projections of the spin along the easy axis. 

The second term in (Q]i, E(S 2 — S 2 ), is the lowest-order anisotropy term that 
can break the uniaxial symmetry of the molecule, and represents a rhombic dis- 
tortion with hard axis x. The Hamiltonian including this term no longer commutes 
with S z , and its eigenstates are now symmetric and antisymmetric superpositions of 
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the I + m), | — m) states. At very low temperatures, such that only the two lowest- 
energy states are thermally occupied, the giant spin can be effectively described as a 
two-level system (that is, a qubit) with eigenstates \S^), \ separated by a tunnel- 
ing splitting 24o [Fig. EJb)] . A giant spin prepared in one of the "classical states" 
\3f ± ) = 2~ 1//2 (|^) ± \s^)), corresponding to the magnetization pointing along ±z, 
would tunnel between | + 3?) and | — 3f) at a frequency h /2Aq in the absence of 
external fields. 

The third term, — gjl^S ■ B, describes the coupling to a magnetic field. If B _L z, 
this represents an extra non-diagonal term which has the effect of rapidly increas- 
ing the tunneling splitting 2Aq(B±) [Fig. [2jc)]. Conversely, B || z has the effect 
of breaking the degeneracy of the "classical" states ±|m). If the longitudinal bias 
^ = gllftS z B z is much larger than 2Aq, the spin is effectively localized on one side 
of the barrier, with vanishing probability to tunnel. 

Finally, the term — La-Ia-A^S represents the coupling to the nuclear spin bath. 
Although in some instances the coupling tensor A may be isotropic (e.g. for the 55 Mn 
nuclei in Mn 4+ ions), this is not true in general. Also, strictly speaking an external 
field acts on the nuclei as well with a term — y^B • 1^ (y is the nuclear gyromagnetic 
ratio), but this can be usually neglected in comparison with the hyperfine coupling. 
If one were to assume that the nuclear spins are static, then the effect of the hyperfine 
field on the giant spin would be the same as that of a static external field, with typical 
strength 10~ 5 — 10~ 3 T, or 0.1 — 10 mK in terms of coupling energy. The spread of 
possible hyperfine couplings is incorporated in a parameter Eq. By inspecting Fig. 
|2jc) for the case B = 0, we see that a transverse component of the hyperfine field 
would have hardly any effect on the nanomagnet's tunneling splitting, whereas a 
longitudinal component gives a bias <^n that is easily several orders of magnitude 
larger than 2Aq. Under these conditions, the giant spin should be frozen for the 
eternity. The case could be made even stronger by accounting also for the dipolar 
interaction between molecules in the crystal, which adds an extra (static) random 
field of order 10~ 3 - 10~ 2 T. So why do we observe tunneling in the experiments? 

The answer comes from the dynamics of the nuclear spins. At the very least, the 
nuclear spin bath will have an intrinsic dynamics due to the mutual dipolar cou- 
plings, which can be described by an additional term of the form Y,k,k'^kVk,k'^-k'- 
This effectively generates a bias on the giant spin that fluctuates on a timescale of 
the order of the nuclear T2 (~ 1 — 10 ms typically). The amplitude of this fluctuation 
can be sufficient to sweep the electron spin levels through the tunneling resonance, 
for that tiny minority of molecules that finds itself having £ ~ at some time. The 
sweeping of hyperfine bias through the resonance yields a Landau-Zener - like in- 
coherent tunneling probability. Once a giant spin has tunnelled, two crucial things 
happen: (i) the distribution of internal dipolar fields in the crystal suddenly changes, 
giving other molecules a chance to have B, ~ 0, etc. . . ; (ii) the hyperfine field on the 
nuclear spins belonging to or surrounding the tunnelled molecule suddenly changes, 
stimulating further dynamics of the nuclear spins. The consequences of (i) are the 
formation of a "hole" in the distribution of dipolar biases 142L with width related 
to the spread of hyperfine bias, and a«v^ l aw f° r the short-term relaxation of the 
magnetization |33|. Here we concentrate on the aspect (ii), namely what happens 
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(a) s (b) (c) 




Fig. 2 (a) Energy levels scheme for the giant spin of Mni2-ac in zero external field, neglecting 
non-diagonal terms in the effective spin Hamiltonian. (b) Including the non-diagonal terms, the 
exact eigenstates become symmetric and antisymmetric superpositions of the "classical" localized 
states, separated by a tunneling splitting 2Aq. (c) Tunneling splitting as a function of external field 
applied perpendicular to the easy axis of anisotropy for Mnl2-ac. The hatched area represents the 
typical range of hyperfine couplings. 



to the nuclear spins as a consequence of a tunneling event. This point is extremely 
interesting and instructive, because it radically deviates from the framework under 
which the dynamics of nuclear spins is commonly analyzed |[Tl [37l . 

Practically all of the theory of nuclear magnetism is based on perturbation the- 
ory, since one typically has a large static magnetic field, B z , applied from the out- 
side, plus some local fluctuations of much smaller amplitude. For instance, one can 
estimate the longitudinal nuclear relaxation rate, , by evaluating the magni- 
tude of the transverse component of the fluctuating field, b±, and assuming that 
the time-correlation of the fluctuations decays exponentially with time constant T, 
(bj_(t\)bj_(t2)) — {b\) exp(— (f2 — t\)/i). Within perturbation theory, this yields: 



2\ 

Tl 1 ^l + G^T 2 ' 



(2) 



with 7 and (On the nuclear gyromagnetic ratio and Larmor frequency, respectively. 
Eq. (0 simply means that the nuclear relaxation rate is proportional to the power 
spectrum of the fluctuating field, taken at the nuclear Larmor frequency. The reader 
shall recognize that this is just the fluctuation-dissipation theorem, a well-known re- 
sult of linear response theory |20|. Expressions like (O are ubiquitous in the NMR 
literature, since they relate an experimentally accessible quantity (the nuclear 7\) 
to the dynamic properties of the environment where the nuclei are immersed (t), 
thereby giving NMR its status of "local probe" for the dynamics of complex sys- 
tems. However, the nuclear spin dynamics in tunneling nanomagnets is one example 
where the use of expressions derived from perturbation and linear response theory 
is incorrect and unjustified. This is particularly true for the nuclei that "belong" to a 
magnetic ion in the molecule core, like 55 Mn or 57 Fe. 

Perturbation theory breaks down here because the hyperfine field produced by 
the electron spins on their nuclei (typically ~ 30 T in Mn and ~ 50 T in Fe) is 
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much larger than any externally applicable field. Therefore, when a spin tunneling 
event occurs, the nuclear Hamiltonian suddenly changes by a large amount, and the 
effect of the "fluctuation" of the electron spin cannot be treated as a perturbation. 
In addition, since the electron spin tunneling events are allowed in the first place by 
nuclear spin fluctuations (at least in the small transverse field limit), it is not even 
possible to treat the nanomagnet as an independent source of "field jumps". The 
dynamics of the nanomagnet and its nuclei must be treated together, since the one 
drives the other and vice versa: again, this is a very uncommon situation in NMR. 

The problem is treated formally in the spin bath theory by writing a master equa- 
tion for the probability P-Mf(t) for the giant spin to remain in the | -ft-) until time t 
OTI . For a nanomagnet alone, that would be trivially: 



When coupled to the nuclear spin bath, the tunneling of S has the effect of "coflip- 
ping" some of the nuclei. Given the details of the spin-bath coupling, one can calcu- 
late what is the "natural" number of nuclei that would flip upon a tunneling event, 
by considering two effects, (i) If the hyperfine field acting on a nuclear spin does not 
exactly invert its direction by 180° upon tunneling (which can happen if the nucleus 
is subject to dipolar fields from different nanomagnets, only one of which tunnels 
at some instant, or if an external transverse field is applied), then the quantization 
axis of the nucleus changes, i.e. the spin precesses around a new axis. Quantum 
mechanically, this is equivalent to a flip of the nuclear spin. This mechanism de- 
pends only on the direction and not on the timescale of the hyperfine field jumps, 
(ii) Even if the hyperfine field changes direction by exactly 180°, the nuclear spin 
may still follow if S flips slowly. The nuclear coflipping probability is then propor- 
tional to (co^/Qq) 2 , where Qq is the "bounce frequency" of the nanomagnet, i.e. the 
frequency of the oscillations of S on the bottom of each potential well. However, 
since typically Qq ~ 10 12 s -1 , this mechanism is usually unimportant. Knowing 
the strength and direction of the individual hyperfine couplings, one can calculate 
the average number, A, of nuclear spins that would coflip with S at each tunneling 
event. Given a certain arbitrary state of the nuclear spin bath at some instant, a tun- 
neling event may require a number M of nuclei to coflip with S in order to conserve 
energy. If M 3> A, such an event is essentially forbidden ("orthogonality blocking" 
Oil). This is accounted for in the theory by writing an effective tunneling matrix 
element, A M ~ A Q X M I 2 /M\, which has its maximum value for M = (and coin- 
cides with the half-tunneling splitting of the nanomagnet alone), and goes quickly 
to zero forM ^> A. Thus Pfft(t) may take any of the possible values Pu(t) obtained 
by replacing Am for Aq in Eq. (3). In addition, since depends on the bias | 

[Eq. ©], which has a hyperfine component <^n that can fluctuate over a range Eq, 
one must also average Pm(? ; % ) over the bias distribution, obtaining: 



P n (t) = 1 




(3) 




(4) 
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<PM(^)> 5 -i = ^exp(-r*0, (5) 

Fm= 7^o' (6) 

where r$ represents the rate for S to tunnel accompanied by the cofiip of M 
nuclear spins. Calling §o the energy scale associated with the flip of X nuclei (the 
"natural" number of coflips), and considering that the highest tunneling probability 
is obtained for M = 0, the leading term for the global tunneling rate becomes: 

N 2 J^exp(-|^ B |/^o), (7) 



where E,g = gfiQS z B tot is the static component of the bias. B tot here is the sum of 
the longitudinal components of an externally applied field and of the dipolar field 
from neighboring nanomagnets. Since the latter is itself time-dependent as soon 
as tunneling events start occurring, this leads to an interesting collective dynamics, 
characterized by a square-root time relaxation ll33ll . Starting from a fully magnetized 
sample, „#(0) = ^# sa t, the short-time behavior of the magnetization is: 



sat 



1 - xA/lkhort, (8) 

VhL~|^=0), (9) 

where £b is the spread of dipolar biases. We see therefore that the microscopic 
properties of the spin bath enter directly in the macroscopic relaxation of an ensem- 
ble of nanomagnets, through the strength of the hyperfine couplings [Eq, in Eq. (0] 
and the cofiipping probability [§o, in Eqs. 0,®]. 

From the opposite perspective, we have seen that to each tunneling event one can 
associate a certain nuclear cofiipping probability. However, for the purpose of com- 
paring theory and experiments, it is convenient to translate this into more common 
NMR language by calculating the nuclear longitudinal (T\) and transverse (T^) re- 
laxation times. The latter is relatively easy to estimate because, unless the tunneling 
rate is made extremely high by applying a strong transverse field, one will generally 
have TT <C -T N , i.e. S remains static on the timescale of the transverse nuclear re- 
laxation. On such a short timescale we therefore recover the applicability of the stan- 
dard perturbative treatments, whereby Ti is determined by the nuclear dipole-dipole 
couplings and can be calculated using the van Vleck method l44l [Tl. Conversely, T\ 
is determined precisely by the tunneling rate of the nanomagnet, which determines 
how often the local field on the nuclear spins changes direction. Quite amusingly, 
this problem was first treated by Abragam (Ref. JT), page 478) as ". . . an example 
that has no physical reality but where the result can be obtained very simply . . . "; 40 
years later, that example has found physical reality in quantum nanomagnets! The 
simple result is that, since the local field changes direction at intervals ~ T, then 

mm 
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rf'-2r N . (io) 

Once again, the point to bear in mind here is that F N is itself strongly dependent 
on the nuclear spin dynamics. Moreover, one should be careful before calling this 
the nuclear spin-/aff/ce relaxation rate, as T^ 1 is normally interpreted. We shall 
come back to this issue in the review of the experimental results, but here we simply 
note that no phonon bath has been introduced so far, in the context of nuclear-spin 
mediated tunneling. A phonon-mediated tunneling rate, F'' 1 , can also be calculated 
HI EH, yielding: 

r*~—*&M— (in 

where Wa is the linewidth of the electron spin states due to intra-well phonon- 
assisted transitions [21]. At low temperatures this tunneling rate is orders of magni- 
tude smaller than the nuclear-driven one. 



3 Experimental results 

From the theoretical treatment of the coupled nanomagnet-nuclear spins system, we 
have obtained several testable predictions. First of all, it is predicted that an initially 
magnetized sample would relax in a non-exponential way, given by Eq. ((HJ. This has 
been verified very early on by magnetization measurement on Fes at very low tem- 
peratures, which show T-independent relaxation below T ~ 360 mK, and indeed, 
a square-root time relaxation function |46| [Fig. |3ja)]. Second, but perhaps more 
important for the present discussion, the nuclear-driven tunneling rate F N should 
depend on the details of the hyperfine couplings, through the spread of hyperfine 
bias, Eq, and the energy scale associated with the "natural" number of nuclear spins 
that coflip with the nanomagnet upon each tunneling event, §q. Both can be changed 
by isotopic substitution of some constituents. Iron-based nanomagnets, like Feg, are 
particularly suited for this type of study because the strength of the spin-bath cou- 
plings can be either increased by replacing 56 Fe (7 = 0) with 57 Fe (/ = 1/2), or 
decreased by replacing 'H (/ = 1 /2) with 2 H (/ = 1 but 6.5 times smaller gyromag- 
netic ratio). When repeating the low-T magnetization relaxation experiments on iso- 
topically substituted samples, it was found that the relaxation rate would increase 
with 57 Fe enrichment, and decrease with 2 H substitution |48| [Fig. Ob)]. Notice 
that in both cases the mass of the molecule is increased, so that no phonon isotope 
effect could explain the observed trend. Finally, by measuring the magnetization de- 
cay while applying a special sequence of longitudinal fields, it was confirmed that 
the fluctuating hyperfine bias creates a "tunneling window" within which the nano- 
magnets can undergo quantum relaxation [46, 48, 42]. The width of this tunneling 
window was found to be dependent on the isotopic composition of the sample. Thus, 
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the effect of nuclear spins on the quantum tunneling of an ensemble of nanomagnets 
has been verified very early on and in a rather uncontroversial way. 

The opposite effect, i.e. the influence of quantum tunneling of the nanomagnets 
on the nuclear spin dynamics, was observed short thereafter for l43l and 57 Fe in 
Fes H, and 55 Mn in Mni2 Il25ll27l [Fig. [2c)]. In each experiment, a ^-independent 
nuclear T\ was found below a certain temperature, and both for Feg and Mn^ this 
was comparable to the temperature at which, for instance, the magnetic hysteresis 
loops also became T-independent P31 l6l IB1. This can be qualitatively understood 
from Eq. ( [Tol l, which obviously implies that when the T-independent F N is the dom- 
inant electron spin fluctuation rate, also T^ 1 should be T -independent. At higher 
temperatures, thermally-assisted transitions in the nanomagnet start to play a role, 
and the situation becomes much more complicated. This has led to some controversy 
in the interpretation of the nuclear relaxation data in the thermally activated regime, 
on which we shall not dwell. Focusing instead on the low-T quantum regime, it's 
worth noting that one could look for the other obvious signature of resonant tunnel- 
ing, i.e. that it should be strongly suppressed by a longitudinal external field. Indeed, 
it has been observed that the nuclear T\ in Mni2 at very low temperatures becomes 
almost two orders of magnitude slower with the application of a small longitudinal 
field [25 . 27), as compared to the zero-field case [Fig.[3ld)]. The interpretation of 
the experiments on Mni2-ac is slightly complicated by the fact that this particular 
compound contains a minority of fast-relaxing molecules ll40l l47l which remain 
dynamic down to the lowest temperatures, whereas the majority species have a neg- 
ligibly small tunneling rate in zero field. On the other hand, the fact that a large 
fraction of the sample can remain fully magnetized for long times (months), has 
allowed to measure the nuclear dynamics in Mni2-ac as a function of the sample 
magnetization, and revealed the effect of dipolar coupling between nuclear spins 
(i.e. the nuclear spin diffusion described by the term £ fc y IkVkk'h' El [34)). Since 
spin diffusion can occur only between nuclei subject to the same local magnetic 
field, a demagnetized sample where half of the spin are "up" and half of the spins 
are "down" should have slower spin diffusion by a factor y2, as compared to a 
fully magnetized sample. This has been indeed verified by measuring the transverse 
relaxation rate 1 l25l l27l . 

So far, we have discussed experiments that essentially confirmed the predictions 
of the theory of the spin bath. Not surprisingly, the situation becomes more puz- 
zling when looking at the thermal properties of quantum nanomagnets and spin 
bath, because the theory does not deal with them. Specific heat is obviously the 
main experimental tool to investigate thermal equilibrium (or lack thereof) in an en- 
semble of nanomagnets. It was found early on that the magnetic specific heat in Fes 
would be unmeasurable at very low temperatures and zero external field, because 
the spin-lattice relaxation time of the electron spins becomes much longer that the 
typical timescale of the experiment (~ 10 3 s at most). However, by applying a large 
transverse magnetic field, the tunneling rate could be made large enough to recover 
the equilibrium magnetic specific heat l23l l24l . Importantly, the case of Mni2-ac 
is different in that one could observe the contribution of the nuclear spins, which 
have a large specific heat at millikelvin temperatures. The magnetic specific heat 
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of Mni2-ac reveals a hyperfine contribution which approaches the full equilibrium 
value when a transverse field is applied, but remains at least partially visible even in 
zero external field. This observation seems to imply that the nuclear spins in Mni2- 
ac find a way to thermalize to the phonon bath even at very low-T, and even without 
large transverse fields. The definitive experimental proof of the thermal equilibrium 
in the nuclear spins of Mni2-ac was then found by directly measuring the 55 Mn nu- 
clear spin temperature by NMR techniques l25l l27l . This is a crucial observation 
because the nuclear spins have essentially no direct link to the phonon bath, there- 
fore their thermal equilibration must proceed via the coupling to the electron spins 
and their subsequent spin-lattice relaxation process. However, we are dealing here 
with a temperature regime where the only electron spin transitions are quantum tun- 
neling ones. The thermalization rate T^ 1 is found to be orders of magnitude faster 
than expected from the known phonon-assisted tunneling rate [Eq. (fTTIl l. and is ac- 
tually very close to the observed nuclear-spin mediated tunneling rate. The precise 
value appears to be a matter of sample size, cooling power and thermal contact to 
the refrigerator, and it's quite plausible that T^ 1 ~ F N in the limit of small sample 
and perfect contact to the thermal bath [25). 

If the nuclear spins are found to be in thermal equilibrium, one may argue that the 
electron spins should be in thermal equilibrium as well. Since the latter are mutually 
coupled by dipolar interactions of the order ^ ~ 0.1 K, at very low-T we may ex- 
pect the ensemble of nanomagnets to undergo a transition to a magnetically ordered 
state, provided the timescale involved in finding the collective ground state is short 
enough. Long-range magnetic ordering in molecular magnets was first found in the 
Mng compound [26 28|, which has negligible anisotropy and therefore maintains a 
very fast electron spin-lattice relaxation time down to the lowest temperatures. The 
electron-spin transitions involved in finding the ordered state are somewhat triv- 
ial (no tunneling), but an additional important result was the confirmation that the 
nuclear T\ measured by NMR precisely coincides with the nuclear spin-Zaff/ce re- 
laxation time as measured by specific heat |28|. The first quantum nanomagnet to 
show long-range magnetic ordering was a special type of Mn4 [8] [Fig.Oe)], char- 
acterized by fast tunneling rate in zero field. In this case, the collective ordered state 
is found through quantum relaxation, and the ferromagnetic phase transition causes 
a peak in the specific heat at the temperature where the nanomagnets undergo long- 
range ordering. Also in this system the nuclear spins were found to be in thermal 
equilibrium. That the nuclear spins should play a role in this process is clear from 
the discussions above, and it has been experimentally verified in Fes, again by play- 
ing with isotopic substitutions. While a "natural" Fes sample falls out of thermal 
equilibrium at low-T in zero field, a sample enriched with 57 Fe shows a much larger 
magnetic specific heat that almost reaches the full equilibrium value [9| [Fig. |3jf)] . 
At the lowest temperatures, the hyperfine contribution to the specific heat is also 
revealed. These results prove once again that the thermalization time, the electron 
spin tunneling rate and the nuclear T\ are closely related and essentially belong to 
the same physical phenomenon, i.e. the inelastic tunneling of the nanomagnet spin, 
driven by the dynamics of the nuclear spin bath, where the phonon bath acts as a 
thermostat for both electrons and nuclei. 
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Fig. 3 Sketches of the qualitative features of some spin-bath related phenomena in quantum nano- 
magnets. (a) Square-root time relaxation of magnetization in Feg |46|, Eq. (I). (b) Isotope effect 
in the short-time magnetic relaxation rate in Feg, Eq. (9)> and r -independent relaxation at low-7 
1481 . (c) Crossover to T-independent nuclear relaxation rate T{~ driven by quantum tunneling in 
Fe 8 EUlll and Mn 12 -ac [23, Eq. 0. (d) Longitudinal field suppression of quantum tunneling in 
Mni2-ac as seen by nuclear Tf l |25 27 1. (e) Long-range magnetic ordering in Mn4 and equilib- 
rium hyperfine specific heat at low-r |8|. (f) Isotope effect in the magnetic specific heat of Feg 



4 Open questions and future directions 

After discussing the experiments involving thermalization of the nuclei and elec- 
tron spins, it should be clear that the missing ingredient in the present description 
of nanomagnet + spin bath is the role of the phonon bath in influencing the tun- 
neling transition in the presence of a fluctuating hyperfine bias. To make the point 
absolutely clear, let us ask the question: "How does a nanomagnet know what is 
its most energetically favorable spin direction at the instant when a tunneling event 
can occur, so that it can participate to a long-range ordered state?" In general, each 
nanomagnet is subject to some bias originating from the dipolar field of its neigh- 
bors, plus the hyperfine bias from the nuclei. Only when the two compensate each 
other, can a tunneling event occur. At that moment, however, the total energy of 
the nanomagnet + nuclei is the same regardless of the direction of the nanomagnet 
spin. The difference between the two possible orientations is that one will have the 
electron spin in a favorable direction with respect to the local dipolar field, but the 
nuclear spins pointing against the hyperfine field, and vice versa in for the other ori- 
entation. It seems that the creation of long-range order in the nanomagnets should 
go at the expenses of the nuclear energy, but we know this is not the case since 
both electron and nuclear spins are found to be in thermal equilibrium down to the 
lowest temperatures, and the equilibrium nuclear specific heat is well visible pre- 
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cisely when long-range magnetic ordering is observed. Thus, the theory needs to be 
improved to include the role of phonons in this process [25 1 . 

One reason to stress this point is that the attention of the community is progres- 
sively being shifted towards coherent tunneling processes in quantum nanomagnets, 
for the purpose of quantum information processing. Then we shall be interested in 
modelling, and ultimately measuring, the decoherence rate of the electron spins un- 
der the most favorable conditions. In particular, it has been shown that when consid- 
ering the effect of nuclear spins and phonons on the decoherence rate, one expects 
a "coherence window" where the coupling to the environment is the least destruc- 
tive (39). Crudely speaking, at low tunneling frequencies the nanomagnet will cou- 
ple strongly to the nuclear spins, because their Larmor frequencies are comparable, 
whereas higher tunneling frequencies would increase the coupling to phonon modes, 
leaving an optimal low-coupling point between these regimes. For dense (i.e., undi- 
luted) and crystalline ensembles of nanomagnets, the dipolar couplings are actually 
more important than the hyperfine ones, and the (strongly T-dependent) optimal 
operation point is determined by the crossing between dipolar and phonon deco- 
herence rates l29l . However, in view of the above discussion on the shortcomings 
of the present description of tunneling in the presence of spin and phonon baths, 
one may wonder to what extent are these predictions accurate. But this is precisely 
the beauty and the strength of the research on quantum nanomagnets: theoretical 
predictions can be tested qualitatively and quantitatively against the experiment, as 
has been done in the past ten years to understand the incoherent tunneling regime. 
It should be possible to single out each contribution to decoherence by virtue of its 
dependence on temperature, field, and isotopic substitution. While the experiments 
needed to demonstrate coherent control of quantum nanomagnets are extremely de- 
manding, the precious information harvested by studying relaxation and dephasing 
in these systems has the potential to push our understanding of nanometer-sized 
quantum systems to unprecedented levels. 
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